Stretching DNA: Role of electrostatic interactions 
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Abstract 

00 : 

'^_>' The effect of electrostatic interactions on the stretching of DNA is investi- 

q ■ gated using a simple worm like chain model. In the limit of small force there 

^ ■ are large conformational fluctuations which are treated using a self-consistent 

C$ , variational approach. For small values of the external force /, we find theo- 

reticlly and by a simple blob picture that the extension scales as fm where 
' To is the Debye screening length. In the limit of large force the electrostatic 

£h ' effects can be accounted for within the semiflexible chain model of DNA by 

assuming that only small excursions from rod-like conformations are possible. 
1-1 , In this regime the extension approaches the contour length as 1/ \J~~f where / 

y—{ \ is the magnitude of the external force. The theory is used to analyze experi- 

ments that have measured the extension of double-stranded DNA subject to 
■ tension at various salt concentrations. The theory reproduces nearly quanti- 

tatively the elastic response of DNA at small and large values of / and for all 
, concentration of the monovalent counterions. The limitations of the theory 

OO ' are also pointed out. 

ON . 

; 

I. INTRODUCTION 

The technical advances in the manipulation of single molecules has enabled the probe of 
mechanical and relaxational measurement in both equilibrium and non-equilibrium con- 
ditions ||. Among the first of such experiments was the investigation of the elastic response 
of A - bacteriophage DNA (A DNA) molecules subject to tension 0. These experiments 
^ and others have offered a window into the behavior of biological molecules on scales ranging 

from several nanometer to few microns. They also provide an opportunity to understand 
the limits of validity of theories based on generic polymer models such as Edwards model or 
simple model of semiflexible chain such as the worm like chain (WLC) model 

The earliest theories describing the elastic response of WLC subject to tension, which 
appear to be relevant to the experiments of Smith et al. , were due to Fixman and Kovac 
H and Crabb and Kovac ||. The precise experiments made possible by nanomanipulation 
of single molecules ]I],|2j has demanded more accurate theories. In these experiments the 
response to a constant force on a magnetic bead attached to A DNA in a solution of varying 
salt concentration is used to probe DNA elasticity. After these pioneering experiments 
were reported several theoretical studies followed 0^0- These papers showed that, when 
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the salt concentration is sufficiently large, then the simple WLC model subject to tension 
quantitatively reproduces the force-extension curves. In particular, the asymptotic approach 
to full extension at large force, was shown to follow the f~ 1 ^ 2 law, where / is the external 
force. 

Since DNA is highly charged one expects ion effects to be extremely important in de- 
termining the elasticity. It has recently been established that the elastic response and 
persistence length of DNA greatly depend not only on the salt concentration but also on 
the valence of counterions |K| . As in the case of simple polyelectrolytes, systematic theo- 
ries of DNA subject to tension is difficult due to the presence of several competing length 
scales. In particular, at intermediate values of the external force, the interplay of several 
forces, namely intrinsic persistence length, electrostatic repulsion, counterion fluctuation, all 
conspire to determine the conformation of DNA. 

Marko and Siggia || included the effects of electrostatic interaction, within a WLC pic- 
ture of DNA by replacing the intrinsic persistence length by a scale dependent effective 
persistence length p| . Such a description implicitly assumes that the intrinsic persistence 



length is large, and hence only small excursions in DNA from rod-like conformations are 
tolerated. The interaction between the charges is assumed to obey the Debye-Hiickel po- 
tential, Vr)H{r) = \b&~ kt jr. The ionic strength of the solution / is related to the screening 
length Tb = k^ 1 through the relation k 2 = AtiIbI, and the Bjerrum length Ib = e 2 /(ekBT) 
represents the strength of interaction. The scale dependent effective persistence length l e Jf 
varies from l Q + Iqsf ( in length scale r ^> r_o ) to l a (r < r^) , where Iqsf = Ib/^k 2 A 2 



12| , |T3|| with A being the mean distance between charges. Marko and Siggia considered the 



limit when the applied force is so large that the angle fluctuation of tangent vector with 
respect to the direction of applied force is small which in consistent with the assumption of 
the OSF theory |T2j|r5| ]. The chain extension z in the presence of tension approaches L like 



J -1 / 2 , in good agreement with experiment at large force. However, when the conformation 
of DNA is coil-like ( at small force ) the theory shows strong deviation from experimental 
results. 

In order to probe the effects of electrostatic interactions at small values of / and I we 
suggest a mean field variational approach. The mean field model of semiflexible chains 



14,15| successfully describes the conformation of the stretching DNA by a constant force 



at high salt concentration, where the electrostatic interactions are negligible ||. In this 
model, the hard constraint on tangential vector u 2 (s) = 1 is replaced by a global constraint 
< u 2 (s) >= 1 so that the average of the magnitude of a tangential vector u(s) is fixed even 
if the magnitude of tangential vectors fluctuate. This theory provides continuous crossover 
formula for the extension from the small force limit to the large force limit. In small force 
limit the extension z increases linear whereas in large force limit the f~ 1 ^ 2 rule is found. 
The theory also reproduces quantitatively experimental the force-extension curves at high 
salt concentration. 

In this paper, we include electrostatic interactions in the mean field model of semiflexible 
chain. This approach provides an adequate theory of DNA elasticity in the small force limit. 
We compare our results with experiments at various salt concentrations. 
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II. SEMIFLEXIBLE CHAIN UNDER TENSION 



We model DNA as a semiflexible worm like chain (WLC) with contour length L. The 
chain can be parameterized by unit tangent vector u(s) = dr(s)/ds, where r(s) represents 
the position vector in three dimensions at the curvilinear position s. The energy cost for 
bending is characterized by the persistence length l p . The tangent vector for WLC chain 
satisfies the local geometric constraint u 2 (s) = 1 for all s. The enforcement of this constraint 
makes the theory for interacting WLC extremely difficult. Recently, it has been shown |9| 
that one can get reliable results for a number of problems involving semiflexible chains by 
replacing the local constraint u 2 (s) = 1 by a global constraint < u 2 (s) >= 1. The resulting 
theory, which in the absence of interactions reduces to the model for semiflexible chains 
proposed by Lagowski et al. |15||, can be systematically derived from a functional integral 



approach. Here we adopt this mean field model to investigate the effects of tension. 

The probability distribution for isolated WLC chain, which preserves the global con- 
straint < u 2 (s) >= 1 is given by 



P Q = exp[-X [ L dsu 2 (S) - r] / L ds(<9u/<9s) 2 - b(u 2 + u 2 )]. (1) 
Jo JO 

The Lagrange multiplier A = 4 (used to enforce the global constraint < u 2 (s) >= 1) 
is (roughly) inversely proportional to the Kuhn length I of the semiflexible chain, and the 
constant rj is related to the persistence length l p . When the free energy of the non-interacting 
chain is optimized with respect to A and we obtain A = 9/ (8/ p ) and 6 = 3/4. These values 
satisfy the constraint < u 2 (s) >= 1. 

The distribution function of the isolated semiflexible chain under tension is given by 

P Q (f) = exp[- [ L ds\(s)u 2 (S) - T] [ L ds(<9u/<9s) 2 + / f (s) • u(s) - b(u 2 + u 2 )] (2) 
J o Jo J 

If the applied force is constant, A(s) has a uniform value for all s ||. It has been shown 
that this assumption gives a self-consistent solution to the stationary condition ||. The 
optimization of the free energy, in the presence of a constant force leads to the modified 
relation between A = 2Z/3 and rj 

The mean square end-to-end distance of semiflexible chain under constant tension can be 
obtained as 

< R 2 > = ( L ( L < ufs') • u(s") > ds'ds" 



J 



2 f 2 



I2lrj. 3 r .. Vf 



/L-/^(l-exp(-^L)) + YL 2 (4) 
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III. EFFECTS OF ELECTROSTATIC INTERACTIONS:INTERPLAY OF 

LENGTH SCALES 



The WLC under the influence of constant tension yields the correct dependence of the 
extension (along the direction of the force) of / at relatively large value of the inverse Debye 
screening length. However, precise experiments on DNA have demonstrated that the simple 
elasticity model is inadequate to take electrostatic effects into account especially at small 
value of the / and at low salt concentrations. Since DNA is highly charged the response of the 
chain to tension clearly depends on the ionic strength and to the valence of the counterions. 
In the presence of the external force there are four important length scales. They are the 
Bjerrum length lg (= e 2 /AirekpT), the Debye screening length k^ 1 = ro, the persistence 
length l p , and the Pincus length £p = ksTJ f [TO. For monovalent ions k 2 = 8irlbP where p 
is the salt (say NaCl) concentration. The persistence length is l p = l + l e where l Q is the bare 
persistence length and l e is the electrostatic contribution. Although it is well accepted that 
for intrinsically stiff chain l e ~ n~ 2 the dependence of l e on k for flexible chains depends on 



a number of factors ||lT Hl3| , P^ H20i ■ 

Since the conformation of DNA subject to tension and salt can change from being rod-like 
to coil-like a proper description requires a detailed understanding of the scale dependence 
of l e . The presence of several length scales and the lack of simple theories for variation of l e 
with k make the development of scaling theories for elasticity of DNA under tension difficult. 

In certain regimes, however, one can devise a scaling type analysis. At low salt concen- 
trations the Coulomb repulsion between various charge residues leads to stiffening of the 
chain. Consider the limit l D < rp, -C L. In this case the chain can be modeled as being 
flexible on scales larger than rp. Inside the electrostatic blob of size rp, monomers are 
stretched due to the Coulomb repulsion. Now consider the case of small force (£p > r D ). 
The mechanical force orients the chain along the tension axis on scales larger than £p. The 
extension is thus given by 

z ~ &(N/N P ) (5) 

where Np is the number of monomers in a Pincus blob, and N = L/a with a being the 
monomer size. The structure within a Pincus blob is imagined to consists of several electro- 
static blobs of size rp. Within each Pincus blob the mechanical forces are not relevant and 
the arrangement of electrostatic blobs is dictated by thermal fluctuations. Thus, 

i P ~ r D (N P a/r D y (6) 

Combining these Eqs. @ and (H), we get 

z/L ~ ti Vv rt- X (7) 

where v is the usual Flory exponent. For a highly charged object, the excluded volume 
interactions are negligible, and hence it is appropriate to take v — 1/2 whence we get at 
small force, 

z/L ~ (f/k B T)r D (8) 
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The blob picture is not appropriate for intrinsically stiff chains under a large external 
force. This region corresponds to £p ~ rp <C l p ~ L. In this case, the orientation of the 
tangent vector is correlated over scales larger than that determined by Coulomb repulsion 
or mechanical energy. Thus segments of chain of length rjj are not independent of one 
another. In this case, as suggested by Marko and Siggia ||, we will show that the applied 
force effectively increases the persistence length of the chain. For values of / > f c — 
ksT / (l + Iosf) the chain adopts a highly stretched conformation. In the limit of large force 
in high salt concentrations Coulomb interactions are unimportant. As discussed previously, 
the conformation here is entirely determined by the competition between entropic elasticity 
and mechanical energy. With this as the background we develop a detailed theory for dealing 
with electrostatic interactions, at the simple Debye Hiickel level, for DNA under tension. 



IV. SELF-CONSISTENT THEORY: SMALL FORCE REGIME 

In order to take electrostatic effects into account we assume, for simplicity, that DNA 
molecule is uniformly charged and each charged segment of the chain interacts via screened 
Coulomb interactions. The Hamiltonian of an isolated DNA molecule consists of non- 
Coulomb part H Q and the electrostatic energy contribution AH. 

H t = H + AH (9) 

The new probability function P[u] including the electrostatic interaction is 

pr 1 p r 1 f ~ AE \ p r i r f L [ L a "a / ^(-^^(s") -r(s / )|) 

P[u] = Po [u]e*p(— ) = P [u]exp[-. ^ J q ds ds ^ - ] (10) 

where uj is equal to Ib/A 2 , A is the distance between the charges. 

When the applied force is smaller than / < ksT/ro, the orientation of the tangential 
vectors are not correlated at scales r > rp, which implies that the mean square average of 
angle fluctuation < 6(s) 2 > is not small. The chain becomes flexible on large scales and 
then adopts "coil like" conformation, although the chain is still stiff on small length scales. 
We use "coil-like" to imply that relatively large excursions from rod-like conformation are 
possible so that the Gaussian approximation, which would be valid when < 9 2 {s) > is small, 
breaks down. 

For this problem the Gaussian approximation employed to treat the long range interac- 
tion to the WLC model is not valid ||. This can be seen by deviations from experiments 
in the extension - force curves at small forces. (See Fig. (6) in ref. ||.) Here we use a self- 
consistent variational theory to describe the effect of small force on the conformations of 
DNA. 

In order to estimate the size of the charged DNA under tension we follow the uniform 



expansion method introduced by Edwards and Singh |2T[. Accordingly, we write 

H t = H x + B (11) 

where 

B = H Q - H l + AH (12) 



5 



Hamiltonian and Hi corresponds to the non-interacting theory in which / is replaced by a 
effective Kuhn length l\ and AH is a perturbation in Hamiltonian. The appropriate value 
of l\ should satisfy (See Eq.(||).) 



< R 2 >= hL - ij^l - e V l ^ L ) + l ^L\ (13) 

For arbitrary choice of l\ , < R 2 > can be rewritten up to first order in B as 

< R 2 >=< R 2 > 1 _ < BR 2 > x + <B >x< R 2 >! (14) 

where <>i indicates average with weight factor exp(—-j^rp). Self-consistency condition re- 
quires that < BR 2 >i=< B >i< R 2 >i so that < R 2 > to first order in B coincides with 
that computed using a reference . We assume that the Kuhn length / = 2/3 A will be replaced 
by li by the coarse graining processes in such a way that new parameter li satisfies Eq.(f|) 
at given force and persistence length l p . Therefore, we have the following self-consistent 
equation: 



(15) 

The right hand side of Eq. (|T5|) is evaluated using 

< Az 2 (k) > AH =< z 2 >-<z 2 > = - ^ (k ; L ;, f) k=o- < z 2 > (16) 

Lr{K, L, I J 

where <> D indicate the average with weight e _Ho//fcsT , and G(r) is the Green function 
associated with the total Hamiltonian 

G(r,L,f) = r (L)=r D[r(s)]exp(^) (17) 

Jr(0)=O K B 1 

and its Fourier transform is 

G(k, L, f) = J d 3 r(s)exp(-ir • k)G(r, L, f). (18) 

If the applied force is constant along the contour then G(k, L, f) = G(k — i/3f, L) [p2]| . We 
can obtain the mean-square average of the end-to-end distance from G(k, L), 

j^(k,L,f) _ ggfXL) , , . 

<Z > ~~ G(k,L,f) |k= °-~ G(k,L) Ik= - i/3f - (19) 

The correlation function in space is obtained by performing functional integral including 
phase factor e _ikr . If we consider the electrostatic interaction as a perturbation, H t = 
H + AH, AH = V DH (r) = l B e~ Kr /r. 

r d 3 k 8 
G(r, L) = / — < e ikr >„-Ei< (AHe ikr ) n > D (20) 
7 (27r) d „ n! 
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We can write the correlation function G{k) ps G Q (k, L)-n Y^ n =o Gi(k, L) 
G Q (k, L) exp(Gi(k, L)) under the Gaussian approximation. 

G {k,L)G 1 {k,L) = <V DH e ikr > Q 

' " ^S" [ ds' I d 3 q-±— < e ^(r(s")-r(s'))+ik.(r(L)- r (0)) > 



(2tt) 3 2J J J V + ^ 2 
k 2 l 

expl (aL - 1 + e~ aL )l x 

6a 



uV ** Us Ids' fdae-( - * , J exp( ^ , J (21) 



16k 2 J J Jo g^i + a/n 2 gili + a/n 2 ' 



with 



a = ( 27^ )2 

- s') = ^-(a(s" -s')-l + e~ a(s "~ s,) ) 
ba 

g J s " - s ') = —(a(s" -s')-l- 2e aL/2 sinhfs" - s')a/2 cosh(L - s" - s')a/2) 
ba 

where the wave vector q is associated with the momentum transfer via electrostatic inter- 
action and 2 ^ fc 2 is the Fourier transform of VD#(r(s") — r(s')). We have introduced the 
dummy parameter a in order to replace three dimensional integral with respect to q to a 
one dimensional integral in w, The mean square average of the end-to-end distance is 

< Az 2 > AH =< z 2 > +(2 g o(k,L)G;(k,L) + c;(kj L))|k= _.^ f (22) 

We solve Eq.(|15D together with Eqs.fl2lD and (p2|) iteratively to find a new coarse grained 
Kuhn length l\, which is related to the persistence length via Eq.(||]). The correction to 
< z 2 > due to the electrostatic interaction is expected to be always positive since electrostatic 
interactions stiffen the chain. Therefore we would expect the effective Kuhn length l\, which 
varies with salt concentration, to be larger than I. 

In order to analyze the experimental measurements at small force using our theory we 
need the parameters L (the contour length), l a (the intrinsic persistence length of DNA), and 
the effective linear charge density I /A. The values of L and l may be obtained by fitting 
the force-extension curve to the data of Smith et al. JIJ at the monovalent salt concentration 
of lOmM NaCl (k^ 1 w 3.2 nm) using WLC model |3J]. For this condition the WLC gives an 
excellent description of the data because the electrostatic interactions are negligible |§||. 
The best fit is obtained with L = 32.7/im and with the intrinsic persistence length l = 53nm. 
The effective charge density 1/A still remains free parameter. In this paper, we choose 
that 1/A in order to find the best fit with experiments. We will discuss the effect of ion 
condensation in the following section. We use the Bjerrum length = 0.7nm in water at 
room temperature with dielectric constant e = 80. 

Before we present the results of the force-extension it is useful to characterize the variation 
of the electrostatic persistence length of DNA with k. Here l e — l p — l Q , l p = 3/i/4 and 
l Q = 53nm. In Fig.(|l]) we plot l e as a function of k. It is clear that there are two distinct 
scaling regimes of behavior. 



7 



l e ~ K 1 K < K x (23) 



K 2 K > K 



x 



where k x « 0.4nm -1 is roughly the crossover value. This behavior is in accord with recent 
theoretical predictions |~7],^3|. The crossover, at least for this DNA, occurs when l e « l Q 
T%| , which is consistent with the condition Iqsf ~ £ - In the inset to Fig. (]]]), we plot the salt 



concentration dependence of the radius of gyration of DNA. In the low salt concentration 
region, k < O.lnm^ 1 , the radius of gyration R g varies like k _1//2 and for high salt concen- 
tration, R g shows little k dependence. The k^ 1 ^ 2 variation, shown in the inset to Fig.(|]) 
implies v — 1/2 for R g . The implication of this calculation, for our purpose, is that in the 
small k regime l e > l Q and as a result the effects of electrostatic interaction dominate at 
scales of l e . Thus, a more detailed theory described here is required to describe the small 
force behavior of DNA. 

In Fig.(0), we plot force-extension curves for DNA molecule at lOmM Na + , ImM Na + , 
O.lmM Na + ion concentration respectively. The large force fits are done by the calculation 
of the angle fluctuations of the tangential vector using L = 32.7/im and l a = 53nm. (See 
next section.) The small force fits are done with self-consistent mean field approximation 
using the same values for the parameters L and l a . 

The choice of A, for which direct measurements are not available, requires explana- 
tion. We choose the effective charge density 1/A as 2.5nm~ 1 , 1.4nm _1 and l.Onm -1 for 
lOmM Na + , ImM Na + , O.lmM Na + ion concentration respectively. These values give the 
best fit to the data. At large salt concentration 10mMNa + , the interaction range (Debye 
radius) rp = 0.3nm is smaller than the size of each base of DNA. We conclude that coun- 
terion condensation is not relevant in this concentration range. If every base pair carries 
a charge — le, the linear charge density will be 2.94nm~ 1 since the size of base pair is ap- 
proximately 0.34nm. The choice of our effective charge density 1/A = 2.5nm~ 1 at lOmM 
NaCl concentration indicates that most of the counterions dissociate from the monomers 
on the backbone. At smaller salt concentration, namely lmMNa + and 0.1mMNa + , the 
interaction range exceeds the Bjerrum length of aqueous solution, l B = 0.7nm. We expect 
counterions are condensed in the vicinity of DNA, which leads to a reduction in the effective 
charge density. According to the Manning condensation theory p^j25|] , the charge density 



larger than one per Bjerrum length leads to counterion condensation. Therefore we choose 
1/A = lAnm^ 1 w 1/Ib- We used smaller value of 1/A for O.lmM concentration. However, 
all the data points from experiment correspond to the rod-like conformation even in the 
small force regime. We expect the self-consistent theory to be valid even at very smaller 
values of force, / < 10~ 2 pN , and the approximate choice of 1/A is then of the order of 1/Ib- 
It is clear from Fig.(fJ) that there are two different regimes in the chain elasticity as 
the applied force increases. They correspond to the "coil-like" and rod-like conformations 
respectively. At a given ionic concentration, we can observe a plateau in the intermediate 
force regime. The validity of self-consistent mean field calculation is limited before the 
onset of plateau, where "coil-like" conformations dominate. Our calculations provide the 
correct estimation of the chain extension when the extension z is much smaller than the 
total contour length L i.e. z/L < 0.5. The mean square average < R 2 > -1 / 2 with zero force 
is proportional to n~ l l 2 and increases as fn~ l as the magnitude of force increases. These 
results agree with experiments and the blob picture (see Eq.(^) in the small force limit. At 
smaller ionic concentrations, the position of the plateau moves to the smaller values of force. 
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In the following section, we discuss the force-extension curves for a semiflexible chain in a 
rod-like conformation. 



V. ELECTROSTATIC EFFECTS ON STRETCHING STIFF DNA: LARGE 

FORCE REGIME 

In this section, we apply the functional integral method to obtain the force-extension 
relation of an intrinsically rigid chain in the limit of large force. This regime was considered 
by Marko and Siggia || who noted that the effects of electrostatic interactions can be 
absorbed into an effective scale dependent persistence length. Here, we provide a derivation 
of this result using a functional integral approach. 

It appears that when a large force is applied to the chain, the segments of the chain 
do not interfere with each other geometrically, because of tension induced stiffening occurs 
on scales from £p = ksT/f. If this assumption is valid the influence of the electrostatic 
interaction can be treated simply by replacing the persistence length l p by an effective 
persistence length V£> = l Q + Iosf- With this replacement the force-extension curve can be 
easily calculated. Since the backbone is intrinsically stiff it follows that the fluctuation in 
angle 9(s), cos 9(s) = u(s) -u(0), is small, we can expand cos(9(s)) ~ 1 + ^9 2 (s) + ^9 3 (s) + ... 
If we consider only fluctuations of the angle, the expansion with respect to 9(s) guarantees 
the constraint of |u(s) 2 | = 1. The Hamiltonian of the system can be written in terms of 9(s) 
as, 

H l p [L <9u(s) 2 u f L .„ f L . Jn , „. D , „ |f(s) 



ocW f^lf ds + ^f ds " f ds'OUyWG-^s" - s') - [ W±cos6(s)ds 
2 Jo os 2 jo J o Jo ksT 

~9 2 {q)Q-\q)dq (24) 
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Q- 1 (q) = G- 1 (q) + l p q 2 + -^-. (25) 



where 9{q) = J ^ /f^ e iqs 9(s)ds, and in the limit of qj k <C 1 



G-\q) = £[(1 + ~^)ln(l + q^-*) - 1] « _ " {q/K y + 0((q/K)% (26) 



Therefore, the effective persistence length can be identified with the coefficient of the 
quadratic term: V^> — l p + Iosf, Iosf = Ib/^^A 2 . 

The generating function Z for the theory given in Eq.(|24|) is 



Z oc / D[9]exp[-- [ L [ L 9{s")G-\s" - s')9(s')ds" ds' + T ^9(s)ds\. (27) 
J 2 Jo Jo Jo ksT 

The ratio of the extension to the total length < z/L > can be obtained from the generating 
function Z 

< Z >- 1 < 62{S) > - 1 19 9 Z (28) 
~L 2 " 1 _ 2 dfdf Z (28) 
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where the mean square value of 9(s) can be found from Q 

1 roo roo 1 roo 

< e 2 (s) >=— dk dk' < 6{k)6{k') > e is{k ~ k } = — / dke 2isk Q{k). (29) 

27T J —oo J —oo 27T J —oo 

In Fig.(@) we show z/L (Eq.^8|) as a function of force. Our results (See Fig.(§)) are in very 
good agreement with experimental results of Smith et al. [Q] for all salt concentrations. In the 
large force regime, we find < z/L >oc —1/y/J which, of course implies, that is regime DNA 
does behave as WLC It is not surprising that the theoretical results for large forces 

start deviating when z/L < 0.4, and dramatically depart from the experimental results and 
when z/L < 0.2 "coil-like" conformations dominate at small forces. This suggests that for 
small force regime a more elaborate theory, such as the one presented in the previous section, 
is required. 

VI. CONCLUSIONS 

In this paper we have considered the effects of electrostatic interactions on the stretching 
of DNA. It is already well established that the simple elastic model at high salt concentration 
gives an excellent description of the response of DNA to tension. Furthermore, it is clear 
that at small values of the concentration of salt and at small force "coil-like" conformation 
become important and an elaborate theory is required. Here we have shown that when the 
applied force is small, the long range interaction in DNA molecule can be properly treated by 
self-consistent variational mean field approximation. Our theory gives excellent agreement 
with experimental results in this regime. The larger the ionic strength, the larger the overlap 
range in the cross over regime from the mean field calculation (coil conformation) to the 
WLC limit (stretched conformation). If the screening length rjj is large compared to the 
total contour length L (td/L ~ 1), then electrostatic interactions dominate at all scales so 
that the chain is effectively stiffened even when the applied force is small. In this case, the 
self-consistent theory is valid only in the limit of very small force. 

One of the limitations of treating the electrostatic interactions, even at the primitive 
Debye-Hiickel level, is that there is no easy way to choose the linear charge density I /A. The 
value of A is essentially controlled by counterion condensation effects. Since the conformation 
of DNA changes upon addition of salt dramatically a proper theory of describing DNA 
elasticity should include fluctuations due to counterion condensation. In the absence of such 
a theory we have used physical arguments to choose a value for A. A more elaborate theory 
that would treat fluctuations on small scales of the value of A as well as larger length scale 
is required to obtain the line density independently. 

It appears that in the limit of small / the self-consistent variational theory may be 
adequate for all values for k. For large value of k and / the effects of electrostatic interaction 
may be treated using such a scale dependent persistence length ||. It would be desirable to 
have a unified theory that can treat both regimes, and hence the case of intermediate values 
of / and k. Such a theory would require using a more elaborate variational Hamiltonian, 
perhaps, similar to the ones used recently to treat the persistence length of polyelectrolyte 
chains [E3J. 

Despite the success of the theory outlined here it is worth pointing out that the WLC 
can only explain the force-extension curves when the counterion is monovalent. It is only 
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in the presence of monovalent ions that the electrostatic persistence length of DNA displays 
the well accepted dependence on k. Baumann et al. have found that multivalent salt 



ions (counterions) have dramatically different effects on persistence length of DNA. By per- 
forming force measurement using laser tweezers they discovered that not only valence but 
also the shape of counterions profoundly affect elasticity of single DNA molecules. This 
observations and other findings by Baumann et al. |T(| clearly suggest that theories and 



simulations that go beyond the simple Debye-Hiickel theory will be required to provide a 
complete description of the response of DNA to external tension. 
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FIGURES 



FIG. 1. The electrostatic contribution to the persistence length l e as a function of k. It is 
clear that l e exhibits k _1 and k~ 2 dependence in low and high salt concentration respectively. 
The inset shows the dependence of the radius of gyration on the salt concentration, which is 
calculated by solving self-consistent Eq.(|l5|). The value of parameters are L = 32.7ram, l a = 53nm, 
1/A = lAQnm' 1 . 

FIG. 2. Comparison between theoretical and experimental force-extension curves for DNA at 
various salt concentrations, (a) 10mMNa + , (h)lmMNa + , (c)0.1mMA r a + The symbols denote 
experimental results, and the dashed lines represent the functional integral estimation for high 
force limit. The force vs. extension curves for small force regime (solid lines) are calculated 
by numerically solving the self-consistent equation. In panel (c) the solid line is calculated with 
A = l.Onm while the dot-dashed line is computed with A = l.lnm. In all cases L = 32.7/xm and 
l = 53nm. 
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